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^ ■ Abstract 

^ . It is believed that some numerical technique must be employed for the determi- 

■ nation of the system parameters of a visual binary or a star with a planet because 

^ ■ the relevant equations are not only highly nonlinear but also transcendental owing 

^ I to the Kepler's equation. Such a common sense, however, is not true; we discover 

I an analytic inversion formula, in which the original orbital parameters are expressed 

^ i as elementary functions of the observable quantities such as the location of four ob- 

^ ! served points and the time interval between these points. The key thing is that we 

. use the time interval but not the time of each observation in order to avoid treating 
O ■ 

' the Kepler's equation. The present formula can be applied even in cases where the 

Qh' observations cover a short arc of the orbit during less than one period. Thus the 

O I formula will be useful in the future astrometric missions such as SIM, GAIA and 

^ ■ JASMINE. 
■ 
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1. Introduction 

The astrometric observation of binaries gives us a lot of invaluable informations on the 
mass of the component stars and the orbital parameters such as the ellipticity (Danby 1988, 
Roy 1988). In the near future, space missions such as SIM ^ , GAIA ^ and JASMINE ^ find a 
number of binaries nearly within lOkpc. A problem in the astrometry arises from a fact that 
we make a measurement of the angular position of the celestial object, which is projected on a 

^ Space Intcrfcromctry Mission (SIM), http://sim.jpl.nasa.gov/ 

^ Global Astrometric Interferometer for Astrophysics (GAIA), http://astro.estec.esa.nl/GAIA/ 

Japan Astrometry Satellite Mission for INfrared Exploration (JASMINE), http://www.jasmine-galaxy.org/ 
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plane perpendicular to the line of sight. For simplicity we call this plane the observed plane. In 
a case of a binary system, the plane of the Keplerian orbit is inclined with respect to the line of 
sight. As a result, the observed ellipse can be considerably different from the original Keplerian 
orbit. To determine the original orbital parameter, we thus need make a kind of data fitting 
(Olevic and Cvetkovic 2004). It rather consumes time and computer resources, especially 
for a huge amount of data which will be available by the near future missions. Obviously 
it would be preferable to use an inversion formula, which enables us to directly determine 
from the observed quantities the original orbital parameters and the inclination. However, 
the conditional equations which connect the observable quantities with the orbital parameters 
are not only highly nonlinear but also transcendental because of the Kepler's equation. It is 
thus believed that an explicit solution is impossible so that some numerical technique must be 
employed to solve the equations for the orbital parameters (Eichhorn and Xu 1990, Catovic 
and Olevic 1992). Such a common sense in this field, however, is not true as shown below. 

The purpose of this letter is to derive the inversion formula: The key point is that we 
use the time interval between the observations but not the time of each observation. First, we 
determine the observed ellipse from the positional data. Next, together with the time interval 
between observations, all of the ellipticity, the orbital period, the major and minor axes of the 
original Keplerian orbit and the inclination are expressed in the measured quantities. 

2. Inversion Formula 

First, we determine the observed ellipse. Next, we discuss how to determine the original 
orbital parameters and the inclination angle. In this letter, we consider only the Keplerian 
motion of the binary by neglecting the motion of the observer and the galactic motion. 

2.1. Observed Ellipse 

We observe an ellipse on the plane perpendicular to the line of sight. The Cartesian 
coordinates on the plane are denoted by {x,y). A general form of the ellipse is 

ax'^ + (3f + 2-fxy + 26x + 2ey = l, (1) 

which is specified by five parameters since the center, the major/minor axes and the orientation 
of the ellipse are arbitrary. We need make at least five observations to determine the parameters. 
The location of each observed point is (xj,|/j) for i = 1, ■ ■ ■ ,5. Then, we find 
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By using the inverse matrix, we can determine (a, /3, 7, 5,5) in terms of {xi,yi). Henceforth, we 
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choose the Cartesian coordinates {x,y) so that the observed eUipse can be reexpressed in the 
standard form as 

- + ^-1 (31 



where a>b. The elhpticity e is ^J~i--WJo? . 

2.2. Time interval 

Before considering a general case, we study a special case to illustrate a problem due to 
the inclination; we assume that an ellipse is inclined around its major axis by the angle i. The 
length of the major and minor axes is denoted by a and 6', respectively. The length of the major 
and minor axes of the observed ellipse becomes a and b'cosi, respectively. Only by using the 
shape of the ellipse, we cannot determine b' and i separately. This degeneracy can be broken if 
we use not only the location of each point but also its time. More rigorously speaking, we use 
the time interval between points as shown below. We adopt four observational data at the time 
of ^4 > ^3 > t2 > ^1- The star is observed at the location of Pj = {xi,yi) = (acosUj,6sinnj) at each 
time ti for i = l,---,4, where a and b have been determined by fixing the observed ellipse in 
the preceding subsection, and Ui denotes the eccentric anomaly. We assume the anti-clockwise 
motion such that Ui > uj for i > j. All we have to do in the case of the clockwise motion is to 
change the signature of Eq. (5) for the area in the following. We define the time interval as 
tij = ti — tj. The number of these additional quantities ^21? ^32 and ^43 is three, which equals the 
number of the unknown parameters of the period, the orientation and angle of the inclination. 

The original Keplerian orbit is assumed to be parameterized by the length of the major 
and minor axes, ax and bx, and the period T. The ellipticity of the orbit is denoted by ex- 
The focus of the original ellipse is projected onto the observed plane at Pg = {xe,ye)- Figure 1 
shows the geometrical configuration. 

One of the important things on the Keplerian orbit is that the area velocity is constant, 
where the area is swept by the line interval between the focus and the object in the Keplerian 
motion (For instance, Goldstein 1980). Henceforth, the object in the Keplerian motion is called 
the star for simplicity. We should note that the projected focus is not necessarily the focus of 
the observed ellipse. Even after the projection, however, the law of the constant area velocity 
still holds, where the area is swept by the line interval between the projected focus and the star. 
The area swept during the time interval t^j is denoted by Sij. The total area of the observed 
ellipse S is irab. The law of the constant area velocity on the observed plane becomes 

- = ^. (4) 
T f- 

2.3. Determining the ellipticity and the orbital period 

By elementary computations, the area Sij is obtained as 
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Fig. 1. A relation between tlie Keplerian ellipse and the observed one. The line of sight is along the 
z-axis. 

Sij = -ab(ui — Uj — — (sinMj — sinuj) + ^(coswj — cosuj)^ . (5) 
It is noteworthy that Sij is hnear in Xe and ye- By using Eq. (4), we obtain 

(6) 
(7) 
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They are reexpressed as 

As - -A, + = 0, (8) 
a b 

Bs - ^B, + y^B, = 0, (9) 
a 



where 



Ai = t2i sinws + ^32 sinwi - ^31 sinM2, (10) 

A2 = t2i COSM3 + ^32 COSMi - t3i COSM2, (11) 

A3 = t2lU3 + ts2Ui - t3iU2, (12) 

Bi = t32SinM4 + t43SinM2 — t42SinM3, (13) 

B2 = ^32 COS U4 + ^43 COS U2 - t42 COS M3 , (14) 

^3 = ^32^4 + ^43^2 -^42M3- (15) 



These equations are solved as 
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A2BS — A3B2 , , 

which determine the location of the projected focus. 

The original major axis is projected onto the observed ellipse at = {xl^Ul) = 
{a cos UL,b sin ul). The ratio of the major axis to the distance between the origin and the 
focus remains same even after the projection. Hence we find 

Pl = — Pe. (18) 

ex 

The positional vector Pl is located on the observed ellipse given by Eq. (3). Thus we obtain 
the ellipticity as 




eK = \l-^ + '^. (19) 
The parameter G [0,27r) is thus given by 

cosml = -i^^ (20) 

sinuL = -p-. (21) 

The location of the projected focus is given by Eqs. (16) and (17) so that we can 
determine the area Sij from Eq. (5). By using 5*21 for instance in Eq. (4), we obtain the orbital 
period as 

T=^. (22) 

'-'21 

2.4- Determining the major axis and the inclination angle 

First, from the above result, we determine the location of the intersection of the observed 
ellipse and the projected minor axis, denoted by P5 = {xs,ys) = {acosus,bsmus). It is not 
necessary to make an observation of the intersection. The major and minor axes divide the 
area of the ellipse in quarter. Even after projecting the ellipse, the projected major and minor 
axes, which are not those of the observed ellipse, divide the observed ellipse in quarter areas. 
This fact implies that us = ul + which does not mean that Pl is perpendicular to P^. 

Up to this point, we have not specified the angle and orientation of the inclination. Let 
the Keplerian ellipse inclined with the angle z in a way that the angle at the origin between 
the periastron and the ascending node is u called the angular distance of the periastron. We 
can assume that the inclination angle is in [0,7r/2), though it takes a value in [0,7r) in the 
standard context of the celestial mechanics. This is because in the present paper we do not 
specify whether the motion of the star in the Keplerian orbit is prograde or retrograde. Hence, 
rigorously speaking, the ascending node, which we have called above for convenience, may be 
the descending node. In short, our inclination angle is between the line of sight and a unit 
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normal vector to the orbital plane, where there exist two unit normal vectors and we choose 
one such as cos? > 0. Only in this paragraph, we adopt another Cartesian coordinates {x',y') 
so that the ascending node is located on the x'-axis. The periastron of the original ellipse is 
projected at = (x^,?/^) = (ai^- cosci;,axsinu;cosz). Similarly an intersection of the ellipse and 
the minor axis is projected at P5 = {x'g,y'g) = {—bx sinuj ,bK cosu cosi) . The component of the 
vector depends on the adopted coordinates. Therefore, it is useful to consider the invariants 
such as |Pi,|, IP5I and |Pl x Ps|. We find 



Pl\ = dxy cos'^uj + sin cucos^i, (23) 



\Ps\ = bxy sin^ iu + cos^ Lucos^i, (24) 

|Pl X P5I = a/^^A'Cosz. (25) 

From Eqs. (23) and (24), we obtain 

C"^ + = al{l + cosH), (26) 

C^-D^ = a^cos2tusin2z. (27) 



Here we used bx = axy^ — defined 

C=|Pl|, (28) 
^=|Q5|, (29) 

by introducing 

Q5 = ^=E=, (30) 

where Q5 denotes a positional vector on a circle which is made by stretching the observed 
ellipse ax/bx times along its minor axis. Equation (25) is rewritten as 

|Pl X Q5I = «/^cosi. (31) 

By using another expression of P^ = (acosUi,6sinUi) and P^ = {— a suiul, b cos ui), C, 
D and \Pl x Q,s\ ^i-re rewritten as 



C=—^xl + yl (32) 
ex 



1 a^yl + b^ 

^-^V l-el ' ^^^^ 

P^xQs\ = ^£=. (34) 
1 



From Eqs. (26) and (31), we obtain 

cos^i — ^cos2 + 1 = 0, (35) 
where we defined 
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By using C'^ + D'^ > 2CD > 2|Pi x Q^], we can show that 

e > 2. (37) 
Equation (35) is solved as 



cos^ = -(e±Ve'-4). (38) 

We can show that ^ + a/^^ — 4 > 2. Since cosi must be no more than unity, we find the only 
solution as 

1 



cos2 = -(e-Ve-4), (39) 

which satisfies cosz G [0,1). Hence we can uniquely determine the inclination angle i for i G 
[0,7r/2), which has been discussed in the subsection 2.4. Here, by substituting Eqs. (32), (33) 
and (34) into Eq. (36), ^ is expressed in terms of a, b, ck, Xe and z/e- 
Next, Eq. (26) determines the length of the major axis as 



aK = \l , , (40) 



1 + cos^ i 

Finally, Eq. (27) determines the orientation of the inclination as 

cos2u; = 2—. (41) 

aj^sin i 

3. Conclusion 

We derive the inversion formula for astrometric observations of binaries. It is summarized 
as follows. First, we fix the observed ellipse by using the location of five points. Second, we 
choose four points and use their locations and the time intervals between these points. The 
projected focus is determined by Eqs. (16) and (17). The ellipticity is given by Eq. (19). The 
orbital period is determined by Eq. (22). The inclination angle is given by Eq. (39). The 
length of the major axis is computed from Eq. (40). Finally, the orientation of the inclination 
is given by Eq. (41). 

Moreover, our result proves that the mapping between a point in the Keplerian motion 
and the observed point (in time and space) is one-to-one if the number of the observation 
is more than four. For instance, a sixth observed point must satisfy our equations with the 
determined values of the parameters. In practice, the observation inevitably associates errors 
so that we must make a kind of fittings for instance by the least square method. Even in such 
a case, our formula would give likely values of parameters so quickly that we could save CPU 
time for fittings. 



We would like to thank N. Gouda, T. Yano and M. Yoshikawa for useful comments. 
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